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Abstract. The connection between the recoupling scheme of four copies of su(l,l), the 
generic superintegrable system on the 3 sphere, and bivariate Racah polynomials is iden¬ 
tified. The Racah polynomials are presented as connection coefficients between eigenfunc¬ 
tions separated in different spherical coordinate systems and equivalently as different ir¬ 
reducible decompositions of the tensor product representations. As a consequence of the 
model, an extension of the quadratic algebra QR(3) is given. It is shown that this alge¬ 
bra closes only with the inclusion of an additional shift operator, beyond the eigenvalue 
operators for the bivariate Racah polynomials, whose polynomial eigenfunctions are deter¬ 
mined. The duality between the variables and the degrees, and hence the bispectrality of 
the polynomials, is interpreted in terms of expansion coefficients of the separated solutions. 
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1 Introduction 


The connection between group theory, special functions and orthogonal polynomials is an area 
that has been of significant interest for many years now and has yielded many beautiful results 
fundamental in theory and in application. In this work, we give a Lie algebraic description of 
Tratnik’s extension of the Racah polynomials [20]. This description relies on the connection 
between these polynomials and representations of the quadratic algebra associated with the 
‘generic’ superintegrable system on the three sphere [13] and its realization in terms of positive 
discrete series representations of su(l, 1) obtained by Genest and Vinet [4]. The system is given 
by the following superintegrable quantum Hamiltonian 


TT _ A a\ a, 2 a,3 a 4 

H - ~ A S 3 + “2+L2+l^ + L2’ 


2 , 2,2 
S 1 + s 2 + S3 


+ S 4 — 1, 


( 1 . 1 ) 


and is referred to a generic since all non-degenerate, second-order superintegrable systems on 
conformally flat pseudo-Euclidean spaces of three dimensions can be obtained through appro¬ 
priate limits from this system [2] . It admits 6 linearly independent constants of the motion and 
hence is considered superintegrable [18]. 

Remarkably, Genest and Vinet [4] have shown that this system can be realized as the tensor 
product of four copies of positive discrete series representations of su(l, 1). They have used this 
representation to identify the 9 j symbols of the algebra and derive identities for the symbols 
using the separated solutions of the Hamiltonian (1.1) in cylindrical coordinates, which are given 
in terms of Jacobi polynomials. The 9 j symbols are given by rational functions multiplied by 
the vacuum coefficients. 

*This paper is a contribution to the Special Issue on Exact Solvability and Symmetry Avatars in honour of 
Luc Vinet. The full collection is available at http://www.emis.de/journals/SIGMA/ESSA2014.html 
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In this paper, we discuss the same tensor product representation but instead focus on bases 
which separate in spherical coordinates and their associated coupling schemes. We show that 
bivariate Racah polynomials can be obtained as expansion coefficients between two spherical 
coordinate systems, or equivalently, two different coupling schemes. The bispectrality of the 
polynomials is obtained naturally from the model. We show also that there is another set of 
commuting difference operators for the bivariate Racah polynomials which fits naturally into 
the scheme. When this difference operator is added to the algebra, it closed to form a quadratic 
algebra on 9 generators, which we will call QR(9) as it is an extension of the Racah algebra QR(3) 
from the univariate case. 

The remainder of the paper is organized as follows. Section 2 contains the necessary back¬ 
ground material on the two-sphere case and the realization of univariate Racah polynomials as 
expansion coefficients between different bases for three copies of su(l, 1). Section 3 contains 
the main results of the paper, namely the connection between the bivariate Racah polynomials 
and the coupling schemes of four copies of su(l, 1). Section 4 contains the results concerning 
the additional difference operator and the algebra QR(9) as well as an interpretation of the 
bispectrality properties of the Racah polynomials in terms of the tensor product representation. 
Section 5 briefly gives some conclusions and future outlook. 


2 Background 

2.1 Positive discrete series representations of su(l, 1) 

Consider the following operator realization of the positive discrete series representation of su(l, 1) 

K ° = 4 + s2 + ^2) ’ K± = 4 (( s ^ ^) 2 “ ^2) ’ ( 2 - 1 ) 

satisfying the su(l, 1) commutation relations 


[K 0 , K±] = ±K ± , [K-,K + \ = 2K 0 (2.2) 

with Casimir operator 

Q = K% - \{K + ,K_} = v{y - 1), a = - 1)(4j/ - 3). 

Here {A, B} = AB + BA. As is well known, see, e.g., [8, 10], the tensor product of multiple 
copies of such representations will itself be a positive discrete series representation of su(l, 1). 
For example, the tensor product of two such representations would have generators 

k$ 2) = kP + k$\ 11 = 0 ,+,-, 

where the upper indices indicate the factor in the tensor product. A subscript is also appended 
to the variable to distinguish each factor in the representation, i.e., 


K. 


(0 



(i) 

and similarly for K± . 

The tensor product can be decomposed into irreducible components. In the case of the 
tensor product of two representations, the irreducible components are indexed by a non-negative 
integer x such that the total Casimir 


= (nf 2) ) 2 - 1 -{k^ 2 \k { } 2) } 
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takes the values ^ 12(^12 — 1) with u \2 = u\ + 1^2 + x, on each irreducible component. In the 
case of two or more tensor products, the decomposition is not canonical. For example for three 
components, V = V (U <8> V^ <8> U®, with algebra generators 

Af 23) = K& + A f) + Af), M = 0, +, 

it is possible to decompose first by coupling and into components of fA 12 ) with each 
irreducible component indexed by the parameter v \2 and then adding the third component. The 
basis in this coupling scheme would be \v 12 , v) with action of operators 

Q (12 Vl2^} = ^12(^12 - 1)|Z'12,*'), ^12 = Ul + V 2 + X, 

QW 12 , v) = v[y - 1)1^12,^), v = Vi + V 2 + V 3 + N. (2.3) 

Here Q is the Casimir operator of the triple tensor product representation of su(l, 1), called the 
total Casimir operator, given by 

Q = (K< 1231 ) 2 - I{4 I23 Ui m) }' (2.4) 

This operator takes the values as in (2.3), indexed by non-negative integer N, on irreducible 
components. Coupling in the other direction, by reducing first <8> into irreducible 
components indexed by v 23 , depending on non-negative integer n, and then into components 
indexed by N determined by the total Casimir operator, gives the following alternate basis 

Q (23 V23, v) = ^23(^23 - 1)1^23, v), 1/23 = ^2 + ^3 + n, 

QW 23 , v) = v{v - 1)1^23, v), V = Vi + V 2 + V 3 + N. 

Note that from its definition (2.4), the total Casimir operator can be expressed as a linear 
combination of the intermediate Casimirs 

Q = Q( 123 ) =Q( 12 ) + Q( 13 ) + Q( 23 )-Q( 1 )-Q( 2 )-Q( 3 ), = ^ - 1). (2.5) 


The extension to 4 copies of su(l, 1) representations is clear. In this case, the total Casimir 
operator has the form 

q = q( 1234 ) = Y Q {ij) ~ 2 Y^ ® (i) - 

Using the operator realization of the positive discrete series representation of su(l, 1) (2.1), 
we can identify the total Casimir operator Q with the Hamiltonian in (1.1) and the intermediate 
Casimirs with the integrals of the motion for the Hamiltonian. Indeed the Hamiltonian H admits 
6 linearly independent integrals of the motion coinciding exactly with the intermediate Casimir 
operators in the discrete series representation (2.1), namely 

_2 


nfe) — _i h fl —s d ) 2 i i u ' JJi i 

^ - A \ Sl ° s j S l° s i) + 4„2 + Ac ,2 + 


a j s i 


CLi + d; — 1 


4 1 4s 2 

The two-sphere version of (1.1) is 
H ( 123 ) = -A S2 + 

and can be expressed as 


03 

“2 + ~2 + “2 > 
s{ s 2 s 2 


4, 2 


2,2,2 1 

S 1 + s 2 + s 3 — 


( 2 . 6 ) 


H (m = 4 q ( 123 ) + £ 

4 

and hence satisfies [i/l 123 ), q(*j)] = 0 . For the three-sphere model, arising from the tensor 
product of four su(l, 1) representations, the Hamiltonian (1.1) is related to the total Casimir via 

H = 4Q (1234 ), [H,Q {ij) }= 0- 
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2.2 Three products of su(l, 1) and the Racah algebra 

The decomposition of three products of su(l, 1) and its connection with the Racah algebra was 
discovered by Genest, Vinet and Zhedanov [6]. The algebra and its connection with Racah 
polynomials will be fundamental for the following sections, so we recall the results here. 
Consider the following interbasis expansion coefficients 

^23,U2 = (^23,^1^12,^), 
and define the following operators: 

k iPv 23 ,v 12 = -^(^3,HQ (12) |z'i2,z') (2.7) 

k 2 P v 2 3,"12 = -^23,^|Q (23 Vl2,^>- (2-8) 

Here, as in the remainder of the paper, the operators are interpreted as acting on the 

vector space V while the operators kj, and later Kj, as acting on the expansion coefficients. 

To determine the operators k\, we see that Q( 12 ) acts on the right, i.e., on the basis 
|z'i 2 ,z'i 23 ,*') giving 

k l P V 23 ,V\ 2 7 ^ 12(^12 1) P U\ 2 ,1723 1 

but Q( 12 ) will also act on the left, i.e., on the (^ 23 ,^) basis, as 

k iPi/ 23 ,ui 2 ^ ] C(i (/h ^'|^'12, v ), (2.9) 

for some expansion coefficients Ch. Similarly, the operator &2 will satisfy 

k 2Pv 23 ,i'i 2 7^23(^23 1 )R/i2,t'23 1 k 2Pu 23 ,v\ 2 77^(1^23, v\ /T, V) • 

In order to determine these expansion coefficients and to hence use the recurrence rela¬ 
tion (2.9) we will identify the algebra generated by k\ and &2 with the quadratic Racah algebra 
QR(3) [9, 21]. Indeed, computing the commutator of k\ and &)2 leads to a third operator kj, 
defined as 

k 3 = [fci, fe 2 ] - 

The equivalent operator in terms of the intermediate Casimir operators is defined as follows 
R= [Q (12) ,Q {23) ], 
so that 

hPu 23 ,v 12 = ^(^23,^|[Q (12) ,Q (23) ]|^12,^> = ^{^23,v\R\vi2,v). 

By direct computation using the su(l, 1) algebra relations (2.2) and the definition of the inter¬ 
mediate Casimir operators (2.1), we see that R can be expressed as the sum over permutations 
of {1, 2, 3} of cubic terms in the su(l, 1) generators, 

R = Y. 2 e ijk K^K^K^. 


( 2 . 10 ) 
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Here e t jk is the sign of the permutation ( ijk ). This formula holds in general; namely for the 
operators R^k defined via 

Rijk = e ijk [Q { - ij \Q^% (2.11) 

the identity (2.10) holds with R replaced by Rijk • In the notation of this section R = R\ 2 z- 
To complete the algebra relations, it remains only to compute the commutators [fa, fa] and 
[fa, fa], which are determined using the algebra relations for the Q 1 ^' 

[R, Q (12) ] = - 2 (Q ^) 2 - 2 {Q (12) , Q (23) } 

+ 2 (Q + Q (1) + QW + Q(3))Q( 12 ) + 2(Q (1 ) - Q&) (Q® - Q) 

[Q^ 23 \r} = -2(Q( 23 )) 2 - 2{Q( 12 ),Q( 23 )} 

+ 2 (Q + Q (1) + QW + Q( 3 ))Q( 23 ) + 2(Q (2 ) - Q&) (qW - Q). 

These relations can be obtained either through the su(l, 1) algebra relations or, more practically, 
using the models (2.1). In terms of the kj operators, the quadratic algebra relations are 

[fa, fa] = k 3 , [fa, fa] = kl + {fa, fa] + dfa + ei, 

[fa, fa] = k\ + {fa, fa] + dfa + e 2 , ( 2 . 12 ) 

with 

■i = l(Q (1, + e ,2) + <3 (3) + Q <123) ). 

e, = + (Q< 3 > - Q , 2 >)(Q m - Q< 123 >), e 2 = ](Q< 2 > - <jM)(Q< 3 > - Q< 123 >). 

From these algebra relations (2.12) it is possible to determine the action of fa as a shift operator 
in the 173 , or equivalently the n index, as well as the action of fa as a shift operators in the v\ 2 , 
or equivalently the x index. Indeed, this direct approach was used in [13]. However, guided 
by hindsight, we see that this algebra is precisely the algebra QR(3) associated with the Racah 
polynomials and we has make use of this identification to determine the action of the opera¬ 
tors fa and fa in terms of the three-term recurrence relation and eigenvalue equations for the 
Racah polynomials. 

Consider now the Racah polynomials defined by 


r n (a,/ 3 , 7 ,< 5 ;x) 


(a + 1 )n(/3 + d + l)n(7 + 1 )n 


x 4+3 


—n, n + a + (3 + 1, — x, x + 7 + 5 + 1 

a + 1,(5 + ^ + 1,7 + 1 


which are polynomials of degree n in \(x) = x{x + 7 + 5 + 1). For consistency with the following 
sections, we shall parametrize the polynomials as 


r n (/3 1 - fa - 1 ,/% - Pi ~ 1 , ~N - 1 ,N + fa\x) 

- ** h; X-N + ^ 

The polynomials satisfy the following eigenvalue equation 


A(x,/3,N)r n (x) = 


-£(7,ftJV)+[*^+) 


r n (x) 


, fa - fa | , ^ 

= K\n, --- I r n (x), 


(2.13) 
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with 


C(x,P,N) = [B(x)(T x - 1 ) + E(x)(T- 1 - 1 )]R n (x), 
, _ {x + Pi- Po){x + Pi)(x + @2 + N)(N - x) 
(2x + ^i)(2x +/3i + 1) 

F( x _ x(x + Po)(N -X- fix + /3 2 )(N + X + Pi) 

(2x + Pi)(2x + Pi — 1 ) 


K,(n, c) = (n + c){n + c — 1). 


The operator C(x,P,N ) is of the form given in [7]. Note that the function n is related to the 


values of the Casimir operators as follows 

^ 12(^12 - 1 ) = k(x, + v 2 ), ^ 23(^23 - 1 ) = n(n, v- 2 + v 3 ). 


The operators associated with multiplication by the transform variable and the eigenvalue 
equation for the Racah polynomials, 




(2.14) 


gives a realization of the algebra QR(3) (2.12) with the following parameters 



Finally, we are in a position to determine the expansion coefficients P U23Vl2 as Racah polynomials 


and identify the operators k\ and k 2 (2.7), (2.8) with kf and kl} (2.14). The identification is 
accomplished by taking 


(2.15) 


Po = 2i/i — 1, Pi = 2vi + 2 v 2 — 1, p 2 = 2z^i + 2v 2 + 2 v$ — 1. 


The identification is only determined up to a possible choice of conjugation, namely 


kf' = G(x)kiG(x) 1 . 


To determine the necessary gauge, it is important to use the requirement that the basis be 
normalized. Therefore, we chose 


Pv 2 3,^12 = G(x,n)Rn(x) 


with 


G(x,n) Pq^i/^2 Pv 23 ,0 


(~ kSf) x {2u 2 ) x {2ui + 2u 2 — 1) X (N + 2lSi + 2u 2 + 2zz 3 — V) x {v\ + V 2 + ^) x 
(2rl)x(— -ZV — 2^3 + l)x( iy i + v 2 — \) X (N + 2v\ + 2v 2 ) x x\ 
n!(2t , 3 ) n (n + 2v 2 + 2z/ 3 — l) n (iV + 2v 2 + 2zz 3 ) n (—IV — 2v\ + l) n 
(~N) n ( 2 v 2 + 2zz 3 ) 2n (2t'2 ) n (N + 2u\ + 2u 2 + 2 z / 3 — l) n 
{2u 2 + 2vz)n(—N — 2v\ — 2 u 2 + l)_v 

(-N — 2iq + 1 )jv(2i , 3 )jv 


X 
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The action of the operator can then be determined as follows 

-TT 23 ^i 2 = k2G(x,n)R n (x) = G(x,n)kfR n (x), 

and so 

/c2 = G(x, n)k^G(x, n)^ 1 . 

Thus finally, we see the action of Q ( 23 ) on the basis \v\ 2 -, v) is given by 

q( 23 ) = G(x, n)A(x, 0, N)G(x, n)~ l = A(x, 0, N), (2.16) 

where T x ]u\2, u) = \v\2 + \,v) and similarly for T~ 1 , where the notation makes use of the fact 
that 1^12 is linear in x. 

3 Four products of su(l, 1) and bivariate Racah polynomials 

In this section we consider four products of su(l, 1) and determine their relation to the bivariate 
polynomials obtained by Tratnik [20]. As above, we obtain the polynomials through interbasis 
expansion coefficients and identify them by comparing the algebra of the operators with the 
algebra generated by the recurrence operators for the polynomials. 

We take as our initial bases the vectors obtained by coupling first with V ^ to obtain 
components and then coupling this with to obtain components and then finally 

with VG) to obtain a complete basis indexed by u\2, ^123 and 12. depending on non-negative 
integers x\, X2 and N, satisfying 

Q (12 V 12 , ^ 123 , Z') = ^12(^12 - 1 ) 1 ^ 12 , Vl23,v), ”12 = 12! +122+ XU 

Q (123 Vl 2 > ^123, v) = ^ 123(^123 - 1)1 ^ 12 , ^123,^), ^123 = ^1 + ^2 + ^3 + ^ 2 , (3.1) 

QWl2, ^123, V) = 12(l2 - 1)| 12 12 , 12123, V), 12 = 12i + 12 2 + 123 + 12^ + N. 

We will build up the action of each of the intermediate Casimir operators Q^ i] ) on this basis 
using the quadratic algebra generated by these 6 operators plus their four linearly independent 
commutator, the Rijk s (2.11), for example 

fli23 = [<3 (12) ,<3 123 '] = [<? (13) ,C3 (12) ] = [0 <23 >,e‘ 13 >]- 

The remaining algebra relations are obtained through direct, though tedious computations 
using the definitions of Rijk and Q' lJ >. They are as follows 

[Rjki, Q (ij) } = {Q m ,Q w } - {Q w , Q ijk) } - 2(Q w - Q U) ) (Q (k) - Q w ) 

+ 2 (qW) + g(*>)gM + 2 (qW + qW)qW 

- 2(g(7) + gW)g( <fc ) - 2 (gW + Q (k) )Q w , (3.2) 

[Rijk,Q (ij) ] = {Q {ij) ,Q {ik) } - {Q (ij) ,Q w } - 2 (g« - gW) (g(« fc ) - Q {k) ). (3.3) 

Notice that for a fixed ijk, the quadratic algebra relations (3.3) can be expressed via Q( yfc ) as 

[Run,Q {ij) ] = —2 - 2 {qM,qW} + 2(gW fc ) + gW + gw + gW)gW) 

- 2(gW - g(j)) (g(v*) _ q(*0) 

and hence replicate the Racah algebra relations on components with fixed values of the Casi¬ 
mir Q&ty by taking k\ = —Q^/2 and /2. This subalgebra will be utilized to obtain 

the action of the operator Q ( 23 ^ on the basis (3.1). 
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3.1 The basis for {Q( 12 \ Q ^ 23 ^} and univariate Racah polynomials 

Let us now consider the basis for V composed of eigenvectors for the commuting operators (ft 2 T 
and (ft 123 \ namely 1 1^23 , z /123 , v) as well as the action of Q 1 ' 23 ' 1 on the basis \v\ 2 , z/ 123 , z/). In order 
to determine that action of (ft 23 ) on this basis, we note that (ft 23 ) also commutes with (ft 123 ) 
and so the action of (ft 23 ) will be just to shift the index u\ 2 . Indeed the action will be the same 
as the product of three representations and so the operator Q I 23 ) will be exactly the eigenvalue 
shift operator of the Racah polynomials given in (2.16) except that the parameter z /123 is no 
longer the total Casimir and is instead given by z /123 = u i + z /2 + z /3 + x 2 - Thus, the action 
of (ft 23 ) on the basis is given by (2.16) with parameters chosen as in (2.15) with the integer N 
is replaced by x' 2 . The expansion coefficients for these two bases are then 

(z/23, z/ 123 , z/|z/i 2 , z/ 123 , ft = (2 z /2 - 1 , 2z/ 3 - 1, ~x 2 -l,x 2 + 2zq + 2v 2 - l;xi), 


with 


WiiCTni = (^23(0), z/ 123 , z/|z/ 12 (®l), Z/ 123 , Z/) {v2z{n\) , Z/123, Z/|l/12(0), Z/ 123 , ft, 
where the following shorthand is used to show the dependence on the integers x\ and rip 
z /23 fti) = ni + z / 2 + z/ 3 , z/i 2 (azi) = xi + z/i + z/ 2 . 

The operators are then 

Q (12) = k(xi,/ 3 i), Q (23) = (a; a;i cr ril ) 1 / 2 A(xi,^,X2)(w.r 1 o- ri i )“ 1/2 = A(xi,ftx 2 ). 


3.2 The basis for {(ft 23 ), (ft 234 )} and bivariate Racah polynomials 

Since we know already [13] that the representation of the algebra generated by the operators is 
associated with the bivariate Racah polynomials of TYatnik, we use this model to build a rep¬ 
resentation of the algebra. In particular, we will identify the operators <ft* J ) and their linear 
combinations in terms of the eigenvalue and recurrence operators found for the Racah polyno¬ 
mials by Geronimo and Uiev [7]. 

The (normalized) bivariate Racah polynomials are defined to be 


R 2 (ni,n 2 -xi,x 2 ;/3-N) = 


s(ft - ft - l,ft - ft - 1, ~X 2 - l,ft + X 2 ] X\) 


{-N) n+m {-N + ft) n+m (ft - ft)n(ft - ft)m 
x r n2 ( 2m +ft - ft - l,ft - ft - l,ni - N - 1,N + m + ft;x 2 - ni). 


(3.4) 


Note that we have dropped the hat notation from [7] since we have no need for the unnormalized 
Racah polynomials. The polynomials (3.4) satisfy the following eigenvalue equations 


Tf R 2 (n; x; ft N) = -nftni + ft - ft - l)i? 2 (n; x; ft IV), 
£fR 2 (n;x;ft JV) = -(m +n 2 )(ni + n 2 + ft - ft - l)R 2 (n;x; /3; N), 


with 

= fto,( T *.-1) + ftl.o)( T r.‘-1). £f = £ c m( T LT,\-i). 

j,k= 0,±1 

£,(1) _ (xi + ft - ft)(xi + ft)(x 2 + xi + ft ) ft 2 - £l) 

(2xi + ft)(2xi + ft + 1) 

£,( 2 ) _ (xi + ft)(xi -I- ft - ft)(x 2 + Xi + ft)(x 2 + Xi + ft + 1)(1V - x 2 )(N + x 2 + ft) 
(1,1 ) (2xi + ft)(2xi + ft + 1)(2x 2 + ft) (2x 2 + ft + 1) 
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r ( 2 ) _ (.T 1 + Pl){x\ + Pi- /3q)(z 2 - Xl)(x 2 + X 2 + (3 2 ) 

(1,0) (2.X1 + /3i)(2xi +/3i + l) 

(2x 2 (x 2 + P 2 ) + 2N(N + p 3 ) + (fo + l)(/3 3 - 1)) 

(2x 2 + /3 2 — 1)(2x 2 + /? 2 + 1) 

^-,( 2 ) _ (2xi(xi + /3i) + (/3o + l)(/3i — 1))(^2 + xi + /3 2 ) 

(0 ’ 1} ~ (2xi + /3i-l)(2xi + /3i + l) 

(x 2 - XI + |S 2 - Pi)(N - x 2 )(iv + x 2 + /S 3 ) 

(2x 2 +/3 2 )(2x 2 +/?2 + 1 ) 

and the other coefficients defined via the inversion operators Ii(f(xi,x 2 )) = /(—x 1 — Pi,x 2 ) 
and / 2 (/(x 1 , x 2 )) = /(xi, -x 2 - # 2 ) as 

For these operators and more, see [7, Appendix A], 

Suppose that the expansion coefficients can be expressed as 


R(n;x;/3;JV) = {u 23 ,v 234 , H^ 12 , ^ 123 , v) = y/u^o r a R 2 {n-, x; /3; AT). 


In terms of the quantum mechanical system (1.1), these coefficients correspond to the ex¬ 
pansion coefficients between eigenfunctions separated in spherical coordinates associated with 
diagonalizing the operators Q^ 123 ^} and the Hamiltonian, related to Q , into one where 

{Q( 23 ), Q( 234 )} are diagonal. The actions of Q 1 ' 12 ' 1 , Q( 23 \ Q (123 ^ on both of the bases \v 42 , ^12.3, v) 
and 1 1^23, ^234, v) can be determined from the results of the previous section. Using their actions, 
we define the operators K\. K 2 , and K 3 that act on the expansion coefficients as 

KiP(n;x;P\N) = -^23, U 234 ,^|Q ( 12 ) |^12, ^123, v) = ^®i, ^ ^ P(n;x;/ 3 ; IV), 

K 2 P(n\ x; / 3 ; N) = -^(u 23 , P234, v\Q { 23 ) Wi2, ^123, (x, /?, N)P(n; x; / 3 ; AT) 

1 /" P 2 ~ Po\ a An 
= - 2 K ( n i>- 2 -j P( ' n; X; ^ 

K 3 P(n;xi-,p; N) = -^23, ^234, ^|< 5 ( 123 ) |^ 12 , ^123, v) = ~k ^x 2 , ^ -P(n; x; /?; IV). 


The shift operator in xi is the same as above (2.13) with IV replaced by x 2 and is given by 

"#2 — Po\ ( P2~ Po 


Ai (x, P, N) = 


-£?(x,/3;IV) + 


— 1 


with the hat indicating conjugation by square root of the ground state ^/uj x a n . 

The action of the operator Q( 234 ) on the expansion coefficients is similarly determined leading 
to the definition of a fourth operator K 4 satisfying 

K 4 P(n;x; P\ N) = -^(i/ 23 , ^ 234 , ^|Q ( 234 Vi 2 , ^ 123 , v) 

= ~^ K ( n i + n 2,^2 + v 3 + u 4 ) P(n; x; /?; N). 

However, the action of Q*- 234 ^, as well as , on the basis |zq 2 , ^ 123 , v), or equivalently as shift 
operators in the variable xi and x 2 , have not yet been determined. Matching the eigenvalues 
of K 4 leads to the identification 

K 4 P(n\ x; P\ N) = — ^Af (x, p, N)P(n-, x; /?; N), 
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where 


A J(x,/3,iV) = 


-£J(x,/3;JV)+ [hxpL. l 


with the remaining (3 given by 

/?3 = 2 vi + 2 z / 2 + 2 ^ + 2 z / 4 — 1 . 

The final operator to be determined is K§, associated with the action of Q( 34 ). Since Q( 34 l 
commutes with Q^ 12 \ its action on the basis will be to shift the variable x 2 . We hypothesize 
that it is of the form 

K 5 P(n;x;/ 3 -N) = -^ 23 , ^ 234 , v\Q^\vi 2 , ^ 123 , "), 

with 

K 5 P(n-x-,/ 3 ;N) = ~^(u x a n ) 1 / 2 Qi(x, ( 3 , N)(ix x a n )~ 1 / 2 P(n;x- ( 3 ; N), 

fMx,/3, AO = -B(x)T X2 - ^(x)T - 1 + S(x) + £(x) + ( ^ 3 + + - ~ ~ • 

The exact form of B(x.) and E(x) are determined from the quadratic algebra relations, and are 
given by 


B(x) = 

£(x) = 


(x 2 + Xi + /? 2 )(x 2 - .Ti + (3 2 - (3\)(N - x 2 )(x 2 + iV + (3 3 ) 
(2x2 + /3 2 + 1)(2x 2 + / 3 2 ) 

(,X2 - Xi)(x 2 + Xi + /3 2 )(A/~ - ^2 + /g 3 - /g 2 )(./V + X2 + (3 2 ) 
(2x2 + @2 ~ 1)(2X2 + P 2 ) 


Thus, we have realized each of the six linearly independent quadratic Casimirs in terms of 
the following operators, which act on the basis |zq 2 , v 123 , v), 


Q il2) = k(xi,Pi), 

(3.5) 

Q (23) = (w x o- n ) 1 / 2 A^(x,/3, A0(w x cr n r 1/2 , 

(3.6) 

Q (123) = k(x 2,132), 

(3.7) 

Q m) = (w x cr n ) 1 / 2 Af(x,^,N)(o; x cr n ) _1 / 2 , 

(3.8) 

Q {1234) = k(N, fa), 

(3.9) 

Q (34) = (w x (T n ) 1 // 2 0i(x,/3,N)(a; x cT n ) _1/2 . 

(3.10) 


Note that the algebra relations do not close without the addition of the shift operator repre¬ 
senting Q^ 34 ). 

4 Some consequences of the model 

In this section, we derive several results concerning the bivariate Racah polynomials which arise 
as a result of the above theory. 
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4.1 The algebra QR(9) 

By construction, the operators (3.5)-(3.10) give simply another basis for the quadratic Casimir 
operators and so the algebra generated by these operators is identical to the algebra defined 
by (3.2), (3.3). As we shall see in this section, this algebra is an extension of the algebra QR(3) 
with 9 linearly independent generators associated with the bivariate Racah polynomials and 
so we call it QR(9). In order to interpret the algebra as an extension of QR(3), we chose the 
following basis 

I<i = -^k(x i,/?i), K 2 = -^Af(x,/3,y), (4.1) 

K 3 = ~ K (x 2 ,/3 2 ), K± = — (x,/3, IV), (4.2) 

AT 5 = -^(x,/3,1V). (4.3) 

Here we have dropped the normalization and weight factors since they have no bearing on the 
algebra relations. Note that whereas the algebra generate by the operators Q^> had 6 linearly 
independent generators, we have restricted to components with a fixed N and so Q, representing 
the total Casimir operator, will be considered as a constant of the algebra and hence gives a linear 
dependence between the 6 generators (2.5). This is analogous to the univariate case and the 
algebra QR(3). 

As we shall see, there will be several copies of the algebra QR(3) in this algebra. The first 
being the one generated by K x and K 2 since these are simply the operators which act on the 
first factor of i? 2 (n; x; (3, N). In order to obtain this algebra, we include the commutator of K± 
and K ‘2 calling it L\. It would be K 3 in QR(3). The copy of the algebra QR(3) is then given by 

[K l ,K 2 ]=L 1 , 

[K 2 M] = Kl + {K u K 2 } + d\Ko + en, [LuKx] = K\ + {K 1 ,K 2 } + d x K x + e 12 , 
di = ^(Q (1) + Q (2) + Q ( 3 ) -2K 3 ), 

en = ~\{Q {3) ~ Q (2) ) (2 K 3 + Q (1) ), e 12 = J (Q (2) - Q (1) ) { 2 K 3 + Q ^). 

Notice that both K\ and K 2 commute with K 3 so d\, en, and ei 2 act as constant for the 
sub-algebra generated by K\, K 2 and L 2 - 
We continue taking commutators and see 

[Ki,K 3 \ = [K 2 ,K 3 } = 0 

but [K ], K 4 ] / 0 and so we call it L 2 ■ The algebra generated by these operators closes again 
to form a copy of QR(3) but only if the operator K 3 is included in the list of generators. The 
algebra relations become 

[K 1 ,K 3 } = L 2 , 

[K a , L 2 \ = K\ + {KnK A } + d 2 K 4 + e 2 i, [L 2 , K\] = K\ + {K U K A } + d 2 K x + e 22 , 
d 2 = l -( 2 K 3 -Q^ -Q^), 

ear = -^(Q-Q (1) )(2A- 5 + Q (2) ) , e 2 2 = ~ Q < " 2 ' I )(2K 5 + Q). 

The operator K 3 is a shift operator in the variable x 2 only so we expect it to generate, 
with K 3 another copy of QR(3). We shall see in the next section that this is due to the fact that 
it is possible to diagonalize K 3 with univariate Racah polynomials in the variable x 2 ■ Defining 

[K 3 ,K 5 } = L 3 , 
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gives 


[K 5 , L :i ] = Kl + {K 3 , K 5 } + d 3 K 3 + e 3 i, [L 3 , K 3 ] = K\ + {K 3 , K 5 } + d 3 K 4 + e 32 , 
d 3 = l -{Q-2K 1 + Q^+Q^), 

e 3 i = (Q (3) - Q (4) ) (Q + 2-Ki), e 32 = ^(Q - Q (4) ) (2Ad + Q (3) ). 

Again, note that Ad commutes with both K 3 and Ad- 

The final copy of QR(3) is generated by AT 2 and R 5 and it is 


L 4 =[K 2 ,K 5 ], 

[K 5 ,L 4 ] = Kl + {K 2 , As} + d 4 K 4 + e 4 i, [La, K 2 \ = K\ + {K 2 , IQ} + d 4 K 2 + e 42 , 
d 4 = -^(Q (2) + Q (3) + Q (4) - 2AT 4 ), 

e 4 ! = - \ (Q (3) - g (4) ) (Q (2) + 2 K 4 ), e 42 = - A (QW - Q (3) ) (Q (4 > + 2 A 4 ). 

As in the previous cases, Ad commutes with both K 2 and A 5 so acts as a constant in this 
subalgebra. 

The operators L\, L 2 , L 3 and L 4 exhaust the set of linearly independent commutators of 
the Aj’s. The only remaining non-zero commutator is 

[K 3 ,K 4 \ = L 4 + L 3 — L 2 — L\. 

The remaining algebra relations defining QR(9) can, as with the previous ones, be directly 
determined by (3.2), (3.3) but we give a set for example. They are 


[^3,^1] 
[Ad, Li] 


[K^Lt] 


0 , 


\ {{Ki,K 2 } + {Ki,K 4 } + {K 2 , IQ} + {A 3 , IQ} - {IQ, KQ - {IQ, K 5 }) 

q ( 4 ) nC 1 ) n( 2 ) 

+ ^-Ad + ^-(Ad + IQ) + Q^(K 3 + IQ) 

+ \ (g (1) g (2) + q w q (4) + q ( 2 ) q (4) ), 

\ {{Ki,IQ} + {AT 3 , Ad} - {Ad, A 2 } - {Ad, Ad} - {K 3 ,K 4 } - {IQ, A 5 }) 

- ^(Ad + Ad) - ^(Ad + Ad) - |a 2 + ^(g (2) + Q^)Q + ^g (2) Q (3) . 


To summarize, as in the univariate case, the algebra generated by the eigenvalue operators 
and multiplication by the variables of the bivariate Racah polynomials can be closed to form 
a quadratic algebra. But unlike the univariate case, an additional operator not arising as a com¬ 
mutator of two basis operators must be adjoined to the algebra. This operator is a shift operator 
in the second variable x 2 . With this additional operator, the algebra closes and has 9 genera¬ 
tors, the operators Ad through Ad and their 4 linearly independent commutators. The algebra 
contains at least 4 subalgebras isomorphic to QR(3). Finally, we emphasize that this algebra is 
isomorphic to the algebra generated by the operators Q^L) whose algebra is much more com¬ 
pactly expressed via (3.2), (3.3). The form of the algebra in terms of the A'* is given here to 
emphasize that it is a natural extension of the algebra QR(3). 
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4.2 Duality in the model 

As another consequence of the model, we describe the duality of the operators mentioned above 
in terms of the expansion coefficients. Recall, that we have 


(^23,^234,^|Q (12 Vl2,^123,^) = K ^1, ^ + 1 
{v 231 ^234, Z/|Q (123) 1^12, I'm,*') = « (x 2 


2 

P2 + I 


(^23, ^234, ^|^12) I'm, V), 

(*'23,*'234,I'|l'12,I'm,I'), 


(^23, ^234, ^|Q (23 Vl 2 , PL23, I') = Af (im, ^234, I'|l'l 2 , I'm, I'), 
(^23, ^234, ^|Q (234 Vl 2 , I'm, I') = Af (l/ 23 , ^234, ^|^ 12 , I'm, I')- 


The operator can of course act on the basis (^ 23 ,^ 234 ) and will give a shift operator of 

0,±1 on each of the values z /23 and U 234 . Similarly, the operator (jA 123 ^ will act on the basis 
^ 23 ,^ 234 } however it will only shift the value z' 2 34 - We can write these as shift operators in m 
and ni + 712 as 

(^23, ^234, ^|Q (12 Vl 2 , ^123, v ) = Af (^23, ^234, ^1^12, I'm, v), 

(^23, ^234, HQ (123 Vl 2 , I'm, v) = Af (l^, ^234, v\vi 2 , V\23-v) ■ 


These two operators hence form a set of commuting shift operators for the same expansion 
coefficients but with a new set of variables x, degrees h and coefficients (3. In the new variables 
the expansion coefficients are 


(^23(£2), ^234(^1), v\vi 2 {n\ + n2),v\2z{n{),v) = yWRi-^n; x; P;N), 
with action of the operators as 

(^23, ^234, HQ (12 Vl 2 , I'm, v) = A|(l/ 2 3, ^234, ^1^12, I'm, I'}, 

(^23, ^234, ^|Q (123 Vl 2 , I'm, v) = Af <1^23, ^234, v\vi 2 , I'm, I'), 

(^23,^234,^|Q (23 Vl2,^123,^) = « f %2, ^ ^ ^ (^23, ^234, v\vi2, ^123, v) 


(^23, ^234, Z'|Q (234) |l'l2,Z'l23,I') =« f^l, ^ l N ) (^23, 


^234, ^1^12, I'm,*') 


The equations necessary to satisfy this duality are then 


K\Xl, 


/?i + 1 


= k h\ + h 2 


k ^ni + n 2i 
as well as 

Af = Af 


/?3 - Po 


= Kl Xi, 


@3 ~ Po 
2 

^1 + 1 


, k\x 2 


, Kin 1 , 


/?2 + l 

2 

p2 - Po 


= Kin i, 


= « x 2 , 


#2 ~ Po 

2 

P 2 T 1' 


Af = Af 


x ^ 1 K ^1V, 

«(aT, 


fli + 1 

2 

/?i + 1 


-£f 


-£3 


K \0, 

K ( 0 , 


P3 - Po 

2 

P2 - Po 


-Cl 


-c% 


(4.4) 

(4.5) 


(4.6) 

(4.7) 
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where are some operators of the form 

7= E D m(. T Wi-h 

j+k= 0,±1 

Solving (4.4), (4.5) as well as the constant terms of (4.6), (4.7) gives exactly the dual variables 
defined in [7], namely 

xi = hi + h 2 + Pz ~ Po + N - 1 , x 2 = hi + P 2 - Po + N - 1 , 

n\ = x 2 + P 2 + N, n 2 = x\ — x 2 + Pi - p 2 , 

/3o = Po, Pi = Po ~ Pz — 21V + 1, 

& = A) - & - 2 N + 1, /3 3 = Po Pi 21V + 1. 

From these identifications, we are able to identify the shift operators as £” = Cj, in agreement 
with [7]. Thus, we see that the duality and therefore the bispectrality of the bivariate Racah 
polynomials arise from symmetries implicit in the model. 

4.3 The action of and another set of bivariate Racah polynomials 

As mentioned above, the operators K 3 and K$ generate a copy of the algebra QR(3) and so 
we anticipate that they will be associated with Racah polynomials as well. Ignoring initially 
the normalization factors, let us consider simply the set of commuting operators {fl 4 (x, P, N), 
Af (x, p, N)}. The operator fl 4 (x, P, N) has the form of an eigenvalue operator for a univariate 
Racah polynomial in the variable x 2 — x\. However, in this new variable the operator Af (x, /?, N) 
would include shifts in this variable of x 2 — x\ —> x 2 — x\ + 2. With this observation in mind, 
we interpret (x, P, N ) as a shift operator in m 2 = x 2 — x±. Indeed, making the identifications 

fhi = xi, rh 2 = x 2 — x \, 

70 = —2A — 2iq — 2v 2 — 217 — 2zz 4 T 1, 

71 = — 2IV - 2v 2 — 2u 3 - 2^4 + 1 , 

72 = -2N - 2v 3 - 2z / 4 + 1, 

73 = — 2N - 2 z / 4 + 1, 

the two commuting operators can be realized as 

(x, P, N)=(n+ (n + -Cf = Af (m, 7 , AO, 

AZ(x,P,N)= (V+2!±^ U + -Cf = Af(m, 7 ,A0. 

A set of simultaneous eigenvectors for these operators are the normalized, Racah polynomials 
A 2 (m;y; 7 ;AT) satisfying the eigenvalues equations 

AfR 2 (m;y;7; N) = k (^y 2 , ^ ^ ^ ^(m; y;7! N ), 

AfR 2 (m;y; 7 ;IV) = n fy x , ^ ^ ^(m; y; 7 ; N). 

To return to x\ and x 2 we transform to the dual variables, or rather as these were set up already 
as the dual variables, we transform back to the variables y \, y 2 and mi, m 2 . In particular, the 
polynomials can be represented as 

R 2 { m; y; f, N ) = ^2(m; y; 7 ; N), 
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with 


yi = X 2 + N + 2ui + 2u 2 + 2u 3 — 1 
mi = y 2 + 72 + N, 


y 2 = x\ + N + 2vi + 2 u 2 - 1 
m 2 = yi-y 2 + z f2 - 72 , 


(4.8) 


and 


70 — — 2N — 2 u\ — 2^ — 2i/ 3 — 2^4 -\- 1 
72 = — 21V - 2u\ - 2u 2 + 1, 


71 = —2N — 2u\ — 2 v 2 — 2^3 + 1 
73 = — 21V — 2v\ + 1. 


The eigenvalue equations for the polynomials in y, defined via (4.8), and m are 



Thus, we see that the operator f2f(x,/3,iV) is another operator that commutes with 
Af(x, /3,1V) and the set of these operators are diagonalized by the bivariate Racah polynomials 
R 2 (m,y, 7 , Af). 

To return to the model. The expansion coefficients from the basis \ v\ 2 , ^ 123 , v) to |^ 34 , 1734 , v) 
can be expressed as 

(^ 34 (mi), ^ 234(^1 + m 2 ), ^|^i 2 (a:i), ^ 123 (^ 2 ), v) = y/u y a m R 2 (mi,m 2 -, yi,y 2 ) 7 > N ), 

which are Racah polynomials in the variables (4.8). They can be represented as a product of 
one Racah polynomial of degree mi depending on both x\ and x 2 and a second one of degree m 2 
depending on x\ only. 

5 Conclusions 

In this article, we have seen how the bivariate Racah polynomials are related to different su(l, 1) 
coupling schemes, in analogy with the univariate case [ 6 ]. The bivariate Racah polynomials are 
seen to be expansion coefficients between sets of eigenfunctions for the Hamiltonian (1.1) in 
different choices of spherical coordinates. This Hamiltonian arises as the Casimir operator of 
the tensor product of four copies of su(l, 1) [4] and its conserved quantities are seen to be 
represented in terms of intermediate Casimir operators for different coupling schemes. 

As consequences of the model, we have shown that the algebra generated by multiplica¬ 
tion by the arguments of the polynomials along with the eigenvalue operators for the bivariate 
Racah polynomials do not close to form a quadratic algebra, as in the univariate case, unless an 
additional operator is included, beyond the operators arising directly as a commutator of the 
generators. This additional operator is a shift operator in the second variable, x 2 and commutes 
with shift operator Af (or equivalently £f). Including this additional operator, the algebra is 
generated by 5 operators plus their 4 linearly independent commutators. This algebra is isomor¬ 
phic to the algebra generated by the intermediate Casimir operators for the tensor product of 
four su(l, 1) representations. However, taking the basis as in (4.1)-(4.3) makes explicit that the 
algebra is an extension of the Racah algebra QR(3) and indeed contains at least four subalgebras 
isomorphic to QR(3). We have also seen how the symmetry of the Hamiltonian and the repre¬ 
sentation in terms of the tensor product allows for a natural definition of the dual variables and 
bispectrality properties of the polynomials. Finally, from the model and the algebra relations we 
have identified an additional operator Ilf that commutes with Af but is a shift operator in the 
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second variable only. This additional operator is necessary for the quadratic algebra to close. 
The eigenfunctions of {Hf, A|} are also bivariate Racah polynomials but of different arguments 
and with a different choice of parameters. 

Let us finish by mentioning several other possible consequences of this model that would be 
interesting to investigate. The first is to utilize more explicitly the separated eigenfunctions 
of the Hamiltonian (1.1) to derive integral formulas and identities for the Racah polynomials, 
analogous to the results of [4] for the 9 j coefficients of the same model. A similar investigation 
has also been completed for the univariate case [17]. A second avenue is to extend the model into 
arbitrary dimensions. The IV-dimensional analog of the quantum Hamiltonian (1.1) is clear and 
it remains superintegrable, exhibiting N(N — l)/2 linearly, though not functionally, independent 
integrals of the motion via the (2.6). Similarly, it is clear how to extend the algebraic 

interpretation in terms of the tensor product of N copies of su(l, 1) representations. The algebra 
relations would be directly analogous to (3.2), (3.3), having been expressed in a symmetric and 
almost dimensionally independent manner. The remaining open question would be whether 
the Tratnik polynomials [7, 20], already defined for arbitrary numbers of variables, are the 
appropriate IV — 1 variable extensions of the Wilson polynomials realizing this algebra. Just as 
in the case of the bivariate polynomials, we expect that additional operators beyond the 2 sets of 
IV — 1 commuting operators defining the recurrence and eigenvalues relations will be necessary 
to for the algebra to close. 

Another interesting avenue of research would be to consider the limits of these Racah poly¬ 
nomials in terms of the contractions of algebra. Again, in the univariate case, it has been 
shown that the limits of the Racah polynomials are generated by contractions of the quadratic 
algebras [14] and furthermore that all such contractions are in fact generated by Lie algebra 
contractions [12]. For the bivariate case, some results have already been obtained. Equivalence 
classes of second-order superintegrable systems in conformally flat space have been identified [1] 
as have the contractions between them [2]. It remains to identify the corresponding limits of the 
polynomials and the action on the algebras. As an example, the multivariate Hahn polynomials 
have been represented in a similar manner [5] corresponding to the singular harmonic oscillator 
in 3D, which can be obtained from the Hamiltonian (1.1) by an appropriate limit. The connection 
between these contraction limits, their quadratic algebras and different coupling schemes should 
provide a fruitful context for classifying multivariate orthogonal polynomials and their limits. 

Finally, we mention that this analysis would also be interesting to extend to the q -case and to 
consider the appropriate algebra for the bivariate q -Racah polynomials [11] and the connections 
with Leonard triples [3], the universal enveloping algebra U q ( sfe) [19] and Hecke algebras [15, 16]. 
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